Spin Hall Effect and Spin Transfer in Disordered Rashba Model 
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Based on numerical study of the Rashba model, we show that the spin Hall conductance remains 
finite in the presence of disorder up to a characteristic length scale, beyond which it vanishes 
exponentially with the system size. We further perform a Laughlin's gauge experiment numerically 
and find that all energy levels cannot cross each other during an adiabatic insertion of the flux in 
accordance with the general level-repulsion rule. It results in zero spin transfer between two edges 
of the sample as each state always evolves back after the insertion of one flux quantum, in contrast 
to the quantum Hall effect. It implies that the intrinsic spin Hall effect vanishes with the turn-on 
of disorder. 
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The two dimensional (2D) electron systems exhibit in- 
teger and fractional quantum Hall effect (IQHE/FQHE) 
in the presence of strong magnetic field 1 . The precise 
connections between the exact quantization of Hall con- 
ductance, the topological property of the wavefunctions, 
and the current carrying edge states in open systems have 
been well established 1 ! 2 ! 3 ! 4 ! 5 ! 6 . As the hallmark of a QHE, 
the quantized plateaus 1 - appear in the presence of dis- 
order where the longitudinal conductance vanishes with 
the opening of a mobility gap, which ensures a robust 
dissipationless transport regime. In an open system, the 
edge states are the chiral current carrying states free from 
backward scattering, which leads to a precise quantized 
charge transfer between edges of the sample as demon- 
strated through Laughlin's gauge argument 3 * 4 ^. 

The recent proposals of the intrinsic and dissipation- 
less spin Hall effect (SHE) in three-dimensional p-doped 
semiconductors^ and 2D electron systems with Rashba 
spin-orbital coupling (SOC) — have stimulated a great 
deal of interest 9 . 10 ! 11 . 12 ! 13 . 14 ! 15 ! 16 ! 17 ! 18 ! 19 . Such intrinsic 
SHE may provide a new way to manipulate electron spins 
in nonmagnetic semiconductors without the application 
of magnetic fields. It was also suggested that the re- 
cent experimentally observed spin polarization or accu- 
mulation 16 in electron systems might be related to this 
effecOLia. 

However, there exits a strong debate regarding the fate 
of SHE in the presence of disorder. A challenge to SHE 
in the Rashba model comes from analytical perturbative 
calculations involving vertex corrections, where SHE is 
found to vanish in the presence of any weak disorder 
On the other hand, numerical calculations based on the 
Kubo formula™ using continuum model in momentum 
space found that spin Hall conductance (SHC) is finite 
for finite-size system with weak disorder, while its proper 
thermodynamic limit is not clear. A consistent picture 
for the SHE is still absent. Especially the topological 
property of the SHE has not been well addressed, while 
all conventional quantum Hall systems are considered to 
have nontrivial topological origin 5 ^, including the Hall 



effect without magnetic field 20-21 . 

In this paper, we present a numerical Kubo formula 
calculation of SHC based on the Rashba lattice model. 
We show that the finite SHC persists in the presence 
of random disorder scattering for finite-size systems up 
to a characteristic length scale, beyond which SHC de- 
creases and vanishes exponentially. The scaling behav- 
ior of the SHC follows a one-parameter scaling law in a 
fashion similar to the longitudinal conductance in a lo- 
calized system 2 ^. Furthermore, we perform numerically 
a Laughlin's gauge experiment 3 ^ by adiabatically insert- 
ing flux to directly probe the experimentally measurable 
spin transfer and accumulation associated with SHE. We 
find that energy levels do not cross each other during the 
adiabatic insertion of the flux once the disorder is turned 
on, in accordance with the level-repulsion rule. It results 
in zero spin transfer between the edges and zero spin ac- 
cumulation at the boundary, as each state always evolves 
back to itself after the insertion of one flux quantum. 
This is in contrast to the quantized charge transfer in the 
IQHE 3,4 -, which is associated with the nontrivial topolog- 
ical property of quantum Hall systems 5 ^. We thus con- 
clude that the topological SHE vanishes in the disordered 
2D Rashba model. An example of a topological SHE in 
2D electron system 21 *^ 3 , is also briefly discussed. 

We start with a tight-binding lattice model of noninter- 
acting electrons with the Rashba SOC. The Hamiltonian 
is given as follows: 24 ' 25 

H = -t c]cj + ^2 Widjci 
(ij) i 

+ Vso ^ cla x c i+ y - i ^2 cl<T y c t+S: + H.c^j (1) 

Here, = (c^, c^) are electron creation operators with 
t as the nearest-neighbor hopping integral and V so the 
Rashba SOC strength, and Wi is an on-site random non- 
magnetic potential uniformly distributed in the interval 
(-W/2,W/2). 
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where the electron wave function, \I/(i) 
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FIG. 1: Spin Hall conductance <j s h (in units of e/4n) as a 
function of Fermi energy Ef (in units of t), for system size 
N = 24 x 24 at W = OAt and V so = 0.5t. The dashed line 
with open triangles is SHC at periodic boundary condition 
averaged over 200 disorder configurations. The solid line with 
filled triangles is SHC averaged over 10 disorder configurations 
and 16 x 16 different boundary phases. 



We numerically study square samples of side L by using 
exact diagonalization method under a twisted boundary 
condition: ^(i+Lx) = e ie "^(i) and ^(i+Ly) = e ie «*(j), 

*t(») 

two spin components. The linear response SHC at zero 
temperature can be calculated by using the Kubo formula 
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FIG. 2: (a) The disorder and boundary phase averaged SHC 
{o~sh) (in units of e/47r) at Ef — — 3.75t as a function of 
disorder strength W (in units of t) for sizes N — 16 x 16, 24 x 
24, 32 x 32 and 48 x 48. (b) < a sH > for W between 0.05-2.0 
t as a function of sample length L — 16, 24, 32, 40, 48. 
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where Ef is the electron Fermi energy, and ^ m repre- 
sents the m-th eigenstate with energy E m . We use the 
standard velocity operator v = j-[H, r] (r is the posi- 
tion operator of electron) and the spin current operator 
jzspm _ 2{ Wx)(Tz j. ) which measures spin current flowing 
along the x direction with the spin polarization along the 
z direction. 

The calculated <j s h at a weak disorder strength W = 
OAt with 200 disorder configuration average, is shown 
in Fig. 1 (dashed line with open triangles) as a function 
of the Fermi energy Ef with the system size N — 24 x 
24 under periodic boundary condition (PBC) with 6 X = 
8 y = 0. Clearly the fluctuations of the SHC are very 
strong much larger than e/47r, in contrast to IQHE where 
the fluctuations in Hall conductance only appear between 
plateaus? 6 with an amplitude smaller than e 2 /h (when 
converting to the units of spin conductance, e 2 /h is equal 
to e/47r). 

The SHC is also very sensitive to the boundary phases 
that we impose at the edges of the finite-size system. For 
example, antiperiodic boundary condition (0 X = 8 y = tt) 
will move the positive and negative peaks in Fig. 1 to 
different Ef's. Remarkably, as we change 8 = (6 x ,6 y ) 
in the 27r x 27t phase space, we find that the abnormal 



positive and negative fluctuations of the SHC at different 
8 can largely cancel each other. Thus (<t s h) averaged 
over 16 x 16 different 0's becomes a smooth function of 
Ef with a value close to 0.3^, which is also plotted in 
Fig. 1 as the line with the filled triangles. Completely 
similar behavior of o s n upon boundary phase average 
has been observed for all W's and V so 's. 

We note that the eigenstates of a finite-size system be- 
come boundary-phase independent only in the thermo- 
dynamic limit. However, when we apply an electric field 
through a time-dependent vector potential, it acts as a 
generalized boundary phase evolving with time. Thus the 
SHC averaged over time is equivalent to the boundary 
phase averaged SHC. Similar boundary phase averaged 
charge Hall conductance gives rise to a topological invari- 
ant Chern numbeB 2 * 2 ^. In Fig. 2(a), we show (<t s h) as a 
function of disorder strength W (in units of t) for system 
sizes N = 16 x 16, 24 x 24, 32 x 32, 48 x 48 with V so = 0M 
at a fixed Fermi energy Ef — — 3.75i near the band bot- 
tom, where the spectrum of the lattice model approaches 
the continuum one. Each data point is averaged over 
random disorder configurations (about 200 for L = 48 
and 500 — 1000 for smaller sizes) and boundary phases 
(16 — 576 different 8 for each disorder configuration). At 
weak W limit, (<j s h) reaches the value of 0.48(^) for all 




FIG. 3: Energy levels of the 2D Rashba model (in units of t) 
with open boundary along y-axis as a function of the twisted 
phase 6 X for a system of V so = O.lt, Ef = — 3.75i and N = 
32 x 32. The solid-line is for W = O.lt and the dotted-line is 
for W = 0. The region indicated by the arrow is enlarged in 
Fig. 4(a). 
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FIG. 4: (a) Enlarged energy levels from Fig. 3 to show the 
avoided level crossing, (b) Robust level crossing in IQHE 
which leads to the adiabatic transfer of edge states. This is 
calculated for the same parameter used in Fig. 3 with an 
additional perpendicular magnetic field. 



sizes, in agreement with the predicted intrinsic value for 
the pure system^. 

The SHC decreases continuously with the increase of 
W as shown in Fig. 2(a). At very small W = 0.05i, the 
SHC is essentially a constant close to the intrinsic value 
with negligible variation with L. However, the lack of a 
flat region in (<j s h) at W — ► limit suggests that the 
effect of disorder is intrinsically important. Very similar 
results are obtained for other Fermi energies and V^ 's. 
Furthermore, we have also seen a monotonic decrease of 
(u s h) with increasing L for W > O.lt. To quantitatively 
understand such a behavior, we plot {<J s h) in a logarith- 
mic scale as a function of L in Fig. 2(b) at different Ws 
with typical error bar shown. In fact all the data can 
be nicely fitted into straight lines, as shown in Fig. 2(b), 
which indicates an exponential decay law 

< (t sH >= c exp (-L/£ s ) (3) 

where the constant cq is insensitive to W (for the whole 
range of VF's we find Co = 0.46 ± 0.04 except for the bot- 
tom curve with W = 2t). Therefore, the scaling behavior 
of SHC can be fitted into the general one-parameter scal- 
ing theory^: a S H = f{L/£ s ). Here, £ s is a spin- related 
characteristic length scale determined by fitting the data 
in Fig. 2(b) into Eq. (3), which approximately follows 
a power-law behavior £ s = 13.8/(W/i) 1 " 5 . It suggests a 
finite length scale for any weak disorder. Thus SHC re- 
mains finite with the sample size up to f s (it can be as 
large as, say, a few hundreds of the lattice constant at 
weak disorder strength W — O.lt), and is extrapolated 
to zero exponentially in the thermodynamic limit. 

To further probe the spin transfer and accumulation 
associated with SHE, we perform numerically a Laugh- 
lin's gauge experiment^ by adiabatically inserting flux 
to a system which is open along the y direction and pe- 
riodic along the x direction. This is the same geometry 



considered by Laughlin for IQHE and the adiabatic in- 
sertion of one flux quantum is equivalent to varying the 
twisted boundary phase X by 27r (note that electric field 
can also be related to the time dependent 9 X ). By di- 
agonalizing the Hamiltonian (1) under open boundary 
condition along the y axis, at 2000 different 9 x 's, the re- 
sulting energies, E m , around Ef = — 3.75t (N = 32 x 32, 
V so = O.lt and W — O.lt) form solid lines shown in 
Fig. 3. Note that the energy spectrum is symmet- 
ric about 6 X = 7r, thus only the half of the spectrum 
with 6 X < 7T is shown. A careful examination of these 
lines reveals that each E m goes up and down, making 
several large angle turns due to backward scatterings. 
These energy levels never cross with each other, except 
for at 9 X = and tt, where two levels become exactly de- 
generate (Kramers degeneracy). Nontrivial topological 
property22*2^ associated with Kramers degeneracy will 
give rise to 2D delocalizationS^. The region pointed by 
an arrow near the left bottom corner of Fig. 3 is enlarged 
in Fig. 4(a), which clearly demonstrates the non-crossing 
feature. This is a general observation for all systems that 
we have checked with different V so 's, Ws and TV's up to 
N = 200 x 200, in accordance with the level-repulsion 
rule of the disordered system. Then if one follows any 
Kramers degenerated pair of states starting from 6 X = 
to 6 — 2-7T, one will always go back to exactly the same 
pair of states. Namely, after the adiabatic insertion of 
one flux quantum, all states evolve exactly back to the 
starting states. Thus it will result in no spin pumping 
between two edges and zero spin accumulation at the 
open boundary. By contrast, in the absence of disorder, 
all energy levels will simply evolve following the apparent 
straight lines (dotted lines in Fig. 3) and cross each other 
with the increase of X . Then after the insertion of one 
flux quantum, each energy level evolves into a new state 
which leads to spin transfer across the sample. But this 
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is the trivial case of level crossing due to the absence of 
disorder scattering. Therefore, the topology of the states 
becomes completely different with the turn on of a weak 
disorder, although the energies themselves only shift by 
very small amounts as shown in Fig. 3. 

For comparison, we show a nontrivial case of level 
crossing in the IQHE for the same Hamiltonian (1) 
in the presence of a perpendicular magnetic field/flux 
(<1>b = 27r/32 per plaquette 2 ^) . As shown in Fig. 4(b), 
the energy levels do show simple crossing. This is due 
to the chiral symmetry of the edge states, which sup- 
presses any backward scattering. Through inserting one 
flux quantum, the two occupied states below Ef are 
pumped onto states above Ef indicated by the arrow in 
Fig. 4(b), which leads to two electrons transferred across 
from one edge to another, corresponding to the IQHE 
with an — 2e 2 /h. We further note the existence of a 
topological SHE in an electron system of 2D graphene 
with Haldane's model 21 - with SOG— , where quantized 
edge transfers are also detected. 

To summarize, we have shown that the SHC persists in 
the presence of random disorder scattering for finite size 
systems up to a characteristic length scale beyond which 
it decreases and vanishes exponentially, in agreement 
with the analytic calculations involving vertex correc- 
tionsii*i£. Thus the length scale £ s should also be identi- 
fied as the spin-relaxation length defined in Refill Fur- 
thermore, by performing numerically a Laughlin's gauge 
experiment^, we have found that all energy levels can- 
not cross each other and thus must result in zero spin 
transfer between the edges and zero spin accumulation 



at the boundary after the insertion of one flux quantum. 
This is in contrast to the quantized charge transfer in the 
QHE 3 4 , which is associated with the nontrivial topologi- 
cal property of quantum Hall systems^SiSi. Therefore we 
establish the vanishing of the topological SHE in the dis- 
ordered 2D Rashba model beyond a linear response the- 
ory with the nonconserved electron spin. The backward 
scattering which reverses electron momentum and spin 
leading to zero spin transfer, is not properly taken into 
account in the linear response theory for SHC as it breaks 
current continuity condition. While finite SHC in four- 
terminal setup is obtained from LB formula^ by mea- 
suring spin current inside the ideal leads, it is related to 
boundary spin current in mesoscopic length scale, which 
cannot be intepreted as bulk SHE due to the nonconser- 
vation of the spin current. 

Note added: After this paper was finished, we be- 
came aware of a numerical work in which the vanishing 
SHC in the thermodynamic limit for an electron Rashba 
model was also concluded, 28 
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